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A generalized relativistic harmonic oscillator for spin 1/2 particles is studied. The Dirac Hamil- 
tonian contains a scalar S and a vector V quadratic potentials in the radial coordinate, as well as 
a tensor potential U linear in r. Setting either or both combinations E = S + V and A — V — S 
to zero, analytical solutions for bound states of the corresponding Dirac equations are found. The 
eigenenergies and wave functions are presented and particular cases are discussed, devoting a special 
attention to the nonrelativistic limit and the case E = 0, for which pseudospin symmetry is exact. 
We also show that the case U — A — is the most natural generalization of the nonrelativistic har- 
monic oscillator. The radial node structure of the Dirac spinor is studied for several combinations 
of harmonic-oscillator potentials, and that study allows us to explain why nuclear intruder levels 
' cannot be described in the framework of the relativistic harmonic oscillator in the pseudospin limit. 
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I. INTRODUCTION 



The harmonic-oscillator potential for relativistic spin 1/2 particles has received considerable attention by many 
groups. The subject is of broad interest since the relativistic harmonic oscillator plays a role in several areas, namely, 
in nuclear and particle physics. In particular, it is the central potential of the nuclear shell model, and it has also 
| been used as the binding (and confining) two-body potential for quarks, with applications in meson and baryon 

■ spectroscopy. 

A special type of harmonic-oscillator potential is achieved by replacing the linear momentum operator p in the 
Dirac equation by p — iflmujr, where (3 is the usual Dirac matrix. This replacement results in a second order 
differential equation for the upper and lower components of the Dirac spinor containing a quadratic potential in the 
radial coordinate r. Besides the quadratic potential, that second order differential equation contains a constant spin- 
orbit term, meaning that, in the nonrelativistic limit, in which only the upper component survives, the degeneracy 
of the energy levels is different from the one found in the non-relativistic harmonic oscillator. As explained in the 
following section, that potential may arise from a Lorentz tensor interaction in spinor space. In Ref. Q] we give a 
list of references where this kind of potential was initially studied. Moshinsky and Szczepaniak Q christened it Dirac 
oscillator and renewed to a great extent the interest in the topic. The Dirac oscillator has been applied to quark 
confinement and supersymmetry [j| and hadron spectroscopy U . Its group symmetries have been studied in Ref. [j| , 
• i-H , its solutions in 2 + 1 dimensions and coupling to magnetic field were worked out in Ref. , and finally its wave packets 

■ and thermodynamical properties in 1 + 1 dimensions were studied in Ref. Q. 

Another possibility to introduce a harmonic potential in the Dirac equation is by mixing vector and scalar harmonic 
potentials with equal magnitude and sign, aiming at obtaining a quadratic potential in the Schrodinger-like second 
order differential equation for each spinor component @ . Kukulin et al. generalized the problem by considering a 
vector-scalar harmonic potential plus the Dirac oscillator. More recently, there has been a wide interest in relativistic 
potentials involving mixtures of vector and scalar potentials with opposite signs. The interest lies on attempts to 
explain the pseudospin symmetry in nuclear physics. Cheng et al. 10], using a Dirac Hamiltonian with scalar S and 
vector V potentials quadratic in space coordinates, found a harmonic-oscillator-like second order equation which can 
be solved analytically for A = V — S — 0, as considered before by Kukulin |9(, and also for S = S + V = 0. Very 
recently, Ginocchio solved the triaxial, axial, and spherical harmonic oscillators for the case A = and applied it to 
the study of antinucleons embedded in nuclei [TT| . The case £ = is particularly relevant in nuclear physics, since it 
is usually pointed out as a necessary condition for occurrence of pseudospin symmetry in nuclei |12L . 

In this work we shall consider a Dirac Hamiltonian that generalizes the previous ones by simultaneously introducing 
radial quadratic potentials for A, E and a linear radial potential for U, the tensor potential defined in the following 
section. We will study the cases with A = and £ = 0, for which analytical bound solutions do exist, including also, 



as particular cases, the Kukulin potential and the original Dirac oscillator. The eigenenergies and eigenfunctions are 
obtained analytically in the general case. Next we shall analyze particular cases and pay a special attention to the 
nonrelativistic limits and to the case U = E = 0, for which pseudospin symmetry is exact and there are still bound 
states, as opposed to what happens with nuclear mean fields. From our analysis we also show that the case U = A = 
is a more natural way to introduce a harmonic oscillator in the Dirac equation than the usual way with U 7^ and 
E = A = 0. The eigenenergies have the usual degeneracy of the nonrelativistic case and the upper component a 
form similar to the nonrelativistic wave function. Actually, by letting the harmonic-oscillator frequency to become 
small compared with the mass, we obtain exactly the eigenvalues and wave functions of the nonrelativistic harmonic 
oscillator. 

We also present the node structure of the radial wave functions, motivated by a study carried on in Ref. |14| for V 
and S radial potentials vanishing as r — > 00. In order to find reasons for the unusual radial node structures in some 
particular cases, we obtain the relations between the radial nodes of the upper and lower components of the Dirac 
spinors by inspection of their analytical forms, and illustrate those relations by plots of several of those components. We 
draw conclusions regarding the impossibility to describe the so-called intruder states by harmonic-oscillator potentials 
with exact pseudospin symmetry. 

This paper is organized as follows. In Sec.^we present the general Dirac equation with scalar and vector potentials 
with harmonic oscillator form, proportional to r 2 , and a tensor potential, linear in r. We then obtain and discuss 
the solutions of this equation for A = and E = 0, respectively. The Dirac oscillator is presented in Sec. IIII Al In 
Sec. IIII Bl and Sec. IIII Cl we look into the particular cases A = and E = in the absence of the tensor potential. We 
analyze the nonrelativistic limits, showing that such a limit does not exist in the case E = 0, up to first order in the 
oscillator frequency divided by the mass. This result is connected to the pseudospin symmetry, as discussed in more 
detail in Sec. lIVI In the same section we also present an analysis of the intruder states and show the impossibility to 
describe them within the harmonic oscillator in the case of exact pseudospin symmetry. Finally, our conclusions are 
summarized in Sec. 

II. GENERALIZED RELATIVISTIC HARMONIC OSCILLATOR 

The time- independent Dirac equation for a spin 1/2 fermion with energy £, in the presence of a potential, reads 

H D V = , (1) 

where the Dirac Hamiltonian is 



H D = a ■ p + /3m + V . (2) 

In this Hamiltonian, m is the fermion mass, p is the momentum operator, and a and (3 are 4x4 matrices which, in 
the usual representation, take the form 

«-(:;)■ >-('<>-')■ (3) 

Here <x is a three-vector whose components are the Pauli matrices, and I stands for the 2x2 identity matrix. The 
matrix potential V in Eq. J2J) may, in general, be written as a linear combination of sixteen linearly independent 
matrices, classified according to their properties under Lorentz transformations: scalar, pseudoscalar, vector, pseu- 
dovector, and tensor. In the following, we exclude the pseudoscalar and pseudovector potentials, only considering 
spherically symmetric potentials, for which the total angular momentum of the system is constant. In spite of being 
called spherically symmetric potentials, they do not necessarily depend only on the radial coordinate r = |r|, because 
of their matrix structure. Indeed, the following potential 

V(r) = V(r)+(3S(r)+i/3a-r U(r) (4) 

is spherically symmetric in this sense, because it commutes with the total angular momentum J = L + S, where L 
and S are the orbital and spin angular momentum operators, respectively. In the last term, r = r/r, and the radial 
functions in Eq. Q are named after the properties their respective terms have under Lorentz transformations: V 
corresponds to the time component of a vector potential, S is a scalar potential, and U is a tensor potential. It is 
important to point out that the Dirac equation for the potential in Eq. (@J is invariant under spatial inversion, hence 
the eigenstates have definite parity. 



It is worth showing where the tensor character of the last term in Eq. comes from. The interaction Lagrangian 
for the tensor coupling of a Dirac spinor ^ with an external field is 



AM 



(5) 



where </> M „ — d^tfiv — d^cf)^, a^ v = (i/2)[y fl , 7"], and 7^ are the Dirac matrices. Considering only the time component 
of the vector potential tf^, <j)°, the corresponding term leads to the following potential in the single-particle (Dirac) 
Hamiltonian: 



2M 



(6) 



For a static radial potential <j>°, the tensor potential U(r) in Eq. is J7(r) = f/(2M) (d<fi°/dr). This term can be 
used to describe a particle with an anomalous magnetic moment. Indeed, the anomalous magnetic interaction has the 
form —ifj,/3a ■ Vy(r), where fx is the anomalous magnetic moment, in units of Bohr magnetons, and (p is the electric 
potential, i.e., the time component of the electromagnetic vector potential, A^. Therefore, the interaction —imuj(3a.-r 
(notice that (3 and a. are now written in reverse order), which defines the "Dirac oscillator" mentioned above, could 
be conceived as an anomalous magnetic interaction generated by an electrostatic harmonic potential. 

In summary, for this type of potentials, the Hamiltonian, the parity operator, and the total angular momentum 
form a complete set of commuting observables. Accordingly, the eigenstates can be classified by the parity (±), by the 
total angular momentum j, and its third component m, quantum numbers. Under these circumstances it is natural 
to use spherical coordinates and the spinor, which is the solution of the Dirac equation ifljl. can be generally written 



f ± (r) 



y (f) 



\ 



y T 



(7) 



(r) J 



In this expression, J^m are the so-called spinor spherical harmonics. They result from the coupling of the two- 
dimensional spinors to the eigenstates of orbital angular momentum and form a complete orthonormal set. The 
orbital angular momentum quantum numbers I and I refer to the upper and lower components, respectively. The 
normalization of the Dirac spinor in Eq. (JJJ implies the following normalization condition on the upper and the lower 
radial functions: 



so that and should be square-integrable functions. 



(l^l a + l/f f | 2 )dr=l. 



(8) 



Using the operator a: • p in the spherical-polar form, Eq. Q can be brought to the form 



la. ■ r 



d 1 

— + - (1+/3K) 
or r 



+ (3[m+S(r)] +V(r) +ifia-fU(r) \t> = £ * 



(9) 



where we have introduced the operator K — — /3(S ■ L+l) whose eigenvalues are k = ± (j + 1/2) . Here £ refers to the 
4x4 matrix whose block diagonal elements are the Pauli matrices, er, and whose off-diagonal blocks are zero. More 
explicitly, the spin-orbit coupling quantum number k is related to the orbital angular momentum quantum number 
by 



-(! + 1) = - (j + 1/2) ,j = l + l/2 (k< 0) 
/ - + (j + 1/2) , 3 =1 — 1/2 (k > 0) 



aligned spin 
unaligned spin. 



(10) 



Actually, the quantum number k completely determines j, I, I and, thus, the parity (—1)': 

3 = \>A-7>, 



Notice that I is given by the same formula as I, but with —k instead of n. Accordingly, the Dirac spinors in Eq. Q 
can be labeled just with k and m, i.e., 



V 



r 

Mr) 



y- Km (r) 



(12) 



Using the property er ■ fy Km — — 3^- Km , the Dirac equation in Eq. Q may be reduced to a set of two coupled first 
order ordinary differential equations for the radial upper and lower components g K and f K , namely, 



d K TT/ N 

— + - + U{r) 
dr r 



g K {r)= [£ + m-A(r)]f K (r) 



(13) 



d_ 

dr 



U(r) 



Mr) 



■[£ -m- T,(r)]g K (r) , 



where we have introduced the "sum" and the "difference" potentials defined by 

E = U + S and A = V-S. 



(14) 



(15) 



Using the expression for f K obtained from Eq. (|13|) and inserting it in Eq. I|14(l we arrive at the following second 
order differential equation for g K : 



d 2 k(k + 1) 
dr 2 r 2 



A' 



d k 

£ + m-A(r) [dr r ^ 



2 K l M + u'{r)-U 2 {r)\g K {r) 



= -[£ -m-£(r)p + m- A(r)]g K (r) , 



(16) 



where the prime means derivative with respect to r. In a similar fashion, i.e., using again both Eqs. I|13fl and <|14[) . a 
second order differential equation for the lower component is obtained: 



dr 2 



k(k — 1) 



£ — m — S(r) 



d k 

U(r) 

dr r 



2k EM. _ u> {r) _ u^ r) U K ( r ) 



= -[£-m-E(r)][£ + m- A(r)]/„(r). 



(17) 



In this paper we shall consider harmonic-oscillator potentials, meaning that potentials S and A are quadratic in r 
and the potential U linear in r. With all these potentials in place, Eqs. (|16Jl and 1)17(1 have to be solved numerically 
because of the quartic potentials in r arising from the product [£ — m — £(r)][£ + m — A(r)] and the terms with 
derivatives of A and S. However, either for A — or for S = the solutions are analytical since, in the first case, 
there are only quadratic potentials in the second order equation for g K [Eq. l(16Jl ]. and, in the second case, the same 
happens in Eq. (|17|) for f K . The case A = U = allows us to get, in the nonrelativistic limit, the spectrum and wave 
functions of the central nonrelativistic harmonic oscillator, i.e, Eq. I|16|) becomes the radial Schrodinger equation for 
the three-dimensional harmonic oscillator. For [7 = 0, both cases E = and A = correspond to SU(2) symmetries 
of the Dirac equation [THl ITfjj . The case E = is related to the pseudospin symmetry in nuclei. 

The aim of this paper is to study such analytical solutions and, in particular, to draw conclusions on the required 
conditions for the pseudospin symmetry to show up. 

Let us first consider 



E(r) = — mid 2 r 2 



A(r) = , U(r) = muj2r , 



(18) 



where ll>2 (a real number) is the frequency related to the tensor potential U and u>i (a non-negative number |28j) 
the frequency related to the "sum" potential E. Actually, this case was already considered by Kukulin et al. [!j for 
UJ2 > 0. The possibility of Dirac bound states with a linear tensor potential, which can be unbounded from below, is 
due to the appearance of U 2 in Eqs. Hlfijl and p7JI. It is noteworthy that the different signs of ui2 give rise to different 
possibilities of signs for the spin-orbit coupling. Equation 1)16(1 for the upper component takes the form 



d 2 k(k+1) 
dr 2 r 2 



—muj 2 (£ + m) + m 2 uj 2 



(2n-l)mLU 2 + (£ 2 -m 2 ) \g K {r) = 



(19) 



It is convenient to introduce the following new variable and parameters: 



m(£ + m) 2 _ 2 j 



nr 2 = a z r\ (20) 



« = a/t^-^ + u?, (2i) 



A = (2k - l)m^ 2 - (£ 2 - m 2 ) ^ ^ 



which allows us to write Eq. (|19|) in the simpler form 



+2^ - - V - A M») = 0, (_>>!! 



dy 2 dy y 

since k(k + 1) = 1(1 + 1). An asymptotic analysis suggests searching for solutions of the type 

g K (y)=Ae-y/ 2 yV +1 V 2 T(y) (24) 

where T(y) is a function yet to be determined and A a normalization constant determined by Eq. JSJ. Inserting this 
expression back into Eq. I(23|) the equation for T(y) reads 



d 2 / 3 \ d If 3 A 
y ^ + { l+ 2- y ky--2{ l+ 2 + -2 



T(y) = 0. (25) 



The solutions of this equation, which guarantee that Vmiy^oc g K (y) = 0, are the generalized Laguerre polynomials of 
degree n, L^y), where 

U l + U±), (26) 



2 V 2 2, 

P = I + \ ■ (27) 

From Eqs. (|26|l and l|2()|l - H22|l . we can get immediately the eigenenergies, which are discrete since n is an integer 
greater than or equal to zero: 



£ 2 - m 2 - (2k — l)mu>2 = ^2rt + I + —J y 2m(E + m)uj 2 + Atti 2 ilj\ . (28) 

If we square Eq. (|28|l the resulting algebraic equation is quartic in E. For particular values of n and k (remember 
that k determines uniquely I) we can solve this equation with respect to E and get the energy of the level with 
quantum numbers (n, k) or, using the standard spectroscopic notation, the level nlj. With the help of the program 
MATHEMATICA we were able to get analytical expressions for those solutions. In general, it can be shown that 
there are at most two real bound state solutions of Eq. (|28|l . one with a positive or zero binding energy and the other 
one with a negative binding energy. As we will see later, it can happen that only positive binding energy solutions 
exist. When both solutions exist, we consider the one with positive binding energy, i.e., for which £ > m. 

Figure n displays the energy for some levels, as a function of u)2, for a fixed u\ and m. From this figure we see 
that for 0J2 = the energy levels with (n = 1 , 1 = 2) are degenerate with the level with (n = 2 , 1 = 0) and have 
no spin-orbit splitting. This is a manifestation of the (2n + I) degeneracy of the levels which can be inferred from 
Eq. (|28|l when 0J2 is set to zero. We will further elaborate on this point in the following section. The spectrum for 
the first seven levels for cl>2 = 2 and oj\ = 1 is shown in Fig. [21 where we see that they are, in general, nondegenerate. 
This feature of the spectrum can be understood by looking at Eq. 1)28(1 ■ since when the tensor potential is present 
(lo2 7^ 0) we have always a K-dependent term (a spin-orbit term) that removes the (2n + I) degeneracy, which is also 
characteristic of the nonrelativistic harmonic-oscillator spectrum. 

It is instructive to get the nonrelativistic limit of Eq. I|28|l . We obtain this limit by letting u!i/m and u)2jvn become 
very small. In this limit, (£ — m)/m — E/m also becomes very small, and £ + m ~ 2m, so that 

E= Un + l + ^j ^lo 2 +oj 2 + U-iL 2 . (29) 
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Figure 1: Energy levels for A = with oj\ = 2 and m = 10 as a function of o^- 



30 



25 - 



^20 



15 



10 



— n = 

-- n= 1 

■-■ n = 2 

■■■ n = 3 



S 1/2 P]/2 P3/2 ^3/2 ^5/2 f 5/2 ^7/2 



Figure 2: Energy spectrum for A = with w\ = 2, ioi — 1 and to = 10. 



The solution of Eq. l|19f) , replacing y by a 2 r 2 , is 

g K (r) = A exp(- 1 a 2 r 2 ) (a 2 r 2 )^' 2 L l ^/ 2 {a 2 r 2 ) . 
For the lower component, the first order differential equation in Eq. (|13|) allows us to write 



£ 



dr 



TTILU2 r 



9k(t). 



(30) 



(31) 



Using the recursion relations for the generalized Laguerre polynomials (see, for example, Ref. |l7j ) we get, for k < 0, 



and, for k > 0, 



exp a r ) (a r ) 

£ + m P V 2 / V ; 



/ muj 2 



ijL'+^ia 2 r 2 )-2L l + Z/2 {a 2 r 2 



(32) 



f K (r) = ^™v(-\a 2 r 2 ) (a 2 r 2 ) l/2 x 



£ + TO 



+ 1+ ~) (i + Li-v v + (n + 1) (i - <- + Y 2 (« 2 



(33) 



Since the generalized Laguerre polynomials of degree n have n distinct zeros |17| we may conclude that g K has n 
nodes, and the expression for f K suggests that it has n nodes for k < and n + 1 for k > 0, 




Next, let us consider the other case 

S(r) = , A(r) = u\ r 2 , U(r) =mLu 2 r, (35) 

where u>2 is again the frequency related to the tensor potential U, and uo\ is the frequency now related to the 
"difference" potential A. We now start from the equation for the lower component, Eq. (|17|l . and introduce the new 
variable and parameters defined by 

V = y 2 =a r , (36) 



- = (2 K +l)m^ 2 -(g 2 -m 2 ) 
a 2 

The resulting equation for / re in these variables is formally the same as Eq. (|23|l . Following the same steps as before, 
to guarantee that f K (r) vanishes when r — > oo, we must have 



2 V 2 2, 

where n is an integer greater than or equal to zero and I is given by Eq. This last quantum number is called 

pseudoorbital angular momentum quantum number in view of its role in the pseudospin symmetry, which will be 
discussed in Sec. II VI Similarly to Eq. (|10|l for k, we may define k in terms of / and s = 1/2 (pseudospin quantum 
number): 

_ = f + =-(j + l/ 2 )> j = 1 + 1/2 (k<0) aligned pseudospin , . 

K ~ \ I =+(j+ 1/2) , j = l- 1/2 (k > 0) unaligned pseudospin. ' ' ' 



In other words, k = —n is a parameter describing the coupling between the pseudo-orbital angular momentum and 
the pseudospin. 

From Eq. (|39(1 we get the following equation for the energy 



£ 2 - 



(2k - l)muj 2 = \2h + I + 



2m(£ — m)uj\ + Am 2 L 



(41) 



Again, there are at most two real solutions of this equation. For positive uj 2 and for most cases when L02 < 0, it can 
be shown that there is always a positive binding energy solution. This is the solution we will take in the following 
examples. Some energy levels for u>\ — 2 and w 2 in a range from to 1 are shown in Fig. For u>2 — there is 
a (2h + I) degeneracy of what we might call a pseudoharmonic oscillator. For a nonvanishing 0J2 this degeneracy 
is removed as we can also see in Fig. |SJ In these figures we use the spectroscopic notation h lj , where states with 
I = 0, 1, 2, . . . are denoted by s ,p , d, . . ., respectively. The ground state for oj\ = 2 and UJ2 = is the state OS1/2, but 
we see from Fig. 0]that it is replaced by the state Opi/2 as the ground state when u)\ = 2 and uj 2 ^ 0.3. 




Figure 4: Energy levels for E = with u>± = 2 and m = 10 as a function of <j)%. 
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Figure 5: Energy spectrum for S = with u)\ = 2, u>2 = lj and to = 10. 



The nonrelativistic limit of the eigenvalue equation in Eq. ()41(l is again reached by decreasing lo\ and U2 such that 
ui/m <C 1 and 1^2/™ <1. If we solve the quartic eigenvalue equation for i?/ro and expand it in powers of wi/m and 
uj2/m we get, in lowest order of these expansion parameters, 



- 3 / 1 

2n + I + - - sign(w 2 ) ( k - - 



(42) 



This equation, valid for S = and u>2 ^ 0, shows that the effect of the harmonic-oscillator A potential is only of second 
order in a nonrelativistic expansion. This important result sheds light on the relativistic nature of the pseudospin 
symmetry to be discussed later in this paper. 

A comparison between Figs. ^ and 21 shows that the effect of the tensor potential (seen by the variation of 00%) , is 
bigger for the case £ = than is for A = 0. As a consequence, we still see in Fig. [5] a quasidegeneracy for some levels, 
in particular, those with the same value of n and that have I = k (aligned spin). This can be explained by the energy 
dependence on uf/m in Eq. (|42J) (for small values of uji/m and u)ijra\ a number smaller than W2, meaning that in 



this case we are near the usual Dirac oscillator [loi = and u>2 ^ 0) , where the states with aligned spin (j 
and the same n are all degenerate, as we will discuss in the following section. 
The lower component radial wave function is 

1 



f K (r) = B exp( - i a 2 r 2 ) (a 2 r^+^L 1 ^ 2 ^ r 2 ) , 



1 + 1/2) 



(43) 



where B is again a normalization constant, determined by Eq. JSJl. As already commented, the degree of the generalized 
Laguerre polynomial determines its number of zeros, so that h gives the number of nodes in f K . 

The upper component for the set of potentials in Eq. Ij35|l is obtained from the first order differential equation in 
Eq. ifTijl. i.e., 



9k{t) = - 



1 



S- 



d 
dr 



muj2 r 



Mr) 



whose explicit form is given, using again the recursion relations of the Laguerre polynomials, for k < 0, by 



9n{r) = 
and, for k > 0, 

9n(r) 



Bh 



exp 



(-ia 2 , 2 )( a 2 r 2 ) 



(Z+2)/2 



muj2 



L 



i+1/2 (a 2 r 2 )-2 4 +3/2 (5 2 r 2 ) 



(44) 



(45) 



Bh 
£ — m 



exp 



h + 1 + 



± ~2 2 

— a r 
2 



muj2 



s2 _2\V2 



(« 2 r 2 ) 



V, 



1-1/2,-2 2\ , /- , ,1/1 . muJ 2 \ T l-l/2,~2 2^ 

(a r ) + (n + 1) ( 1 + -— - ) L fl+ { (a r ) 



(46) 



These functions are analogous to the g K and f K functions of the previous case, so from the above discussion we can 
conclude that the number of nodes of the radial functions for £ = is such that (remember that k = — k) 



1 



"9 ■■ 



K < 
K > 0. 



(47) 



This is illustrated in the plots of g K and f K in Fig. Adopting the definition of n as the number of nodes of the 
lower component when ll>2 = and £ = 0, we see that we can convert from the notation n lj to n lj for the same state 
by setting n — * n = n, I — * I = I — 1 when k > and n — * h = n — 1, I — ► I = I + 1 when k < 0. 

We have noticed, however, that the radial function node structure depends on the ratio u>i/u>2- If ^2 becomes 
significantly larger than ui±, in the case of A = the structure of the nodes is the one of the normal Dirac oscillator, 
the same as in Eq. l(4T|l . as it will be discussed in the following section. 



III. PARTICULAR CASES 



Now we will present some results for the three particular cases u)\ = 0, 0J2 = (A = and S = 0). 



A. Case uj\ — 



As mentioned before, the expression "Dirac oscillator" applies to this case when u>2 > [2j. Here we also consider 
the possibility that UJ2 be negative. From Eq. (|19|l when u)\ — 0, it is straightforward to conclude that the equation 
for the upper radial wave function of the Dirac oscillator is a Schrodinger-like equation with a harmonic-oscillator 
potential and a spin-orbit coupling term. 




The eigenenergies are readily obtained taking the limit lo\ — > in Eq. 1)28)1 . which yields 



£ 2 = m 2 + TOju^l 



2( N + - ) +sgn(w 2 )(2K-l) 



(48) 



where the quantum number N = 2n + I was introduced. There are, in general (when the right-hand side of the 
equation is positive) two symmetric real solutions for this equation. Here we take the positive energy solution. Note 
that Eq. I)48|) results from the equation for the upper component, g K . We can get the nonrelativistic energy by setting 
loi = in Eq. 129)1 . Therefore, in this limit we obtain for the energy E the quantum mechanical result |w2|(^ + 3/2) 
plus the spin-orbit contribution \(2n — 1)^2- This term removes the nonrelativistic degeneracy related to levels with 
the same N in such a way that changing the sign of w 2 leads to a order reversal of the levels with the same (n, I ) and 
j = I ± 1/2 quantum numbers. For u>2 > 0, depending on whether the quantum number k is positive or negative, the 
energy is given by [see Eq. fTT)|l ] 



£ 2 = to 2 + 4 muj 2 \n + I + 



1 



K > 0, 



(49) 



£- 



4 mu)2n, 



k < 0. 



(50) 



Thus, for a given quantum number n, all states with aligned spin j = I + 1/2 (/t < 0) and arbitrary orbital angular 
momentum are degenerate. This degeneracy, as well as the (n + I) degeneracy for n > 0, is illustrated in Fig. in 
which the spectrum for L02 = 1 is presented. On the other hand, the eigenenergies of the Dirac oscillator Hamiltonian 
may also be determined from Eq. (|41|l . which was obtained from the second order equation for the lower component 
f K , very similar to Eq. (|19|l . Again taking the limit uj\ — > in Eq. (|41|l one obtains 



m\u)2\ 



+ 5 ) +sgn(w 2 )(2 K + :i : 



(51) 



where N = 2h + I is a new quantum number. The term 2k + 1 comes from the spin-orbit interaction for the lower 
component. Depending on whether the quantum number k is positive or negative, the energy, for 0J2 > 0, is given by 



£ 2 = 



+ 4raw2 [ h + I + 



k > 0, 



(52) 



£ 2 = m 2 + Amu2 (h + 1), 



K < 0. 



(53) 



Of course, the eingenvalue equations in Eqs. 1)521) and (|53() should lead to the same energy values as Eqs. Q49H and 
l|5Up. Remembering that, for k < 0, one has 1 = 1 + 1, and, for k > 0, I = I — 1, from the comparison of those two sets 
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Figure 7: Energy spectra for the normal Dirac oscillator (uii = 0) with LO2 = 1 and ?ti = 10. 



of equations one concludes that 

{n„ — 1 , k < 
(54) 
n g , k > 0, 

since the quantum numbers n and n are equal to the number of nodes of the upper and lower radial functions, denoted 
by n g and n/, respectively. Notice that n g > 1 for k < 0. This can also be inferred from the analytical form of the 
radial upper and lower functions for the Dirac oscillator. These components are readily obtained from Eqs. (|30|l . i|32[) . 
and J22J), noticing that, when ui\ = 0, the square of the parameter a, defined in Eq. (|22|1 . is equal to mL02 (0J2 > 0), 
giving 

g K (r) = A exp(- l - a 2 r 2 ) (a 2 r 2 )^ 1 )/ 2 ^ +1 / 2 (a 2 r 2 ) (55) 

for the upper component, and 

Mr) = exp(-i a 2 , 2 )( a 2 r 2 ) ( ' +2)/2 e 3 / 2 (a 2 r 2 ) (« < 0) (56) 



■ Ur) = It^ exp (" 5 fl2 r2 ) ( fl2 r2 )' /2 { n + l + l) L ^ V2(a2 r2) (K > 0) (57) 

for the lower component [we have used the relation L"" 1 (x) — L"(x) = — L"_ 1 (x) [l7|]. 

As we have seen before, the degree of these Laguerre polynomials is identified with the principal quantum number 
of the upper and lower radial functions, and thus with their number of nodes, so that the relations in Eq. Q54JI follow 
immediately 

In this case (0J2 > 0) Eq. (|54() implies a peculiar node structure for the ground state Si/ 2 (k = —1). This state 
has a node (n g = 1) in the upper component, as we can see in Fig. El whereas the nonrelativistic ground state wave 
function is nodeless. 

If 0J2 < 0, the infinite degeneracy referred to above occurs for k > 0, i.e., for nonaligned spins. In addition, the 
relation between the values n g and nf is modified: rij = n g + 1 for n > and rif = n g for k < 0. It is worth noticing 
that the structure of the radial nodes for cj 2 positive (negative) is exactly the same for the case E = (A = 0) 
discussed in the preceding section. This same node structure could again be inferred from the form of the radial wave 
functions when L02 < 0, which can be obtained once more from Eqs. H^Ufl . (|32fl . and (|33H . setting a 2 ~ —mui2- We 
note that now the ground state Si/ 2 has zero nodes in the upper component, and therefore, in the nonrelativistic 
limit, when the lower component disappears, the wave function has the same node structure as the nonrelativistic 
harmonic-oscillator ground state wave function. Furthermore, from Eq. I|51|l we see that, for L02 < 0, the sign of the 
spin-orbit interaction is inverted, such that states with nonaligned spin (k > 0) have now lower energy than states 




with the same I and aligned spin. This is the situation in atomic physics, where, for example, the states P3/2 have 
higher energy than the states Pi/2- The opposite happens in nuclear physics, where states with aligned spin are lower 
in energy than non-aligned spin states, which in the present case happens when 0J2 > 0. 




(a) r (b) r (c) r 



Figure 9: Radial wave functions for the normal Dirac oscillator of the states (a) 2si/ 2 , (b) W3/2, and (c) 1^5/2 with 
the same parameters as in Fig. 

We present in Fig. El the radial wave functions for the states 2si/2, 1^3/2, and 1^5/2- From these plots we see 
that the upper components for the spin-orbit partners (same n,l) are equal [see Fig.|^Jb) and Fig.[5Jc)], since both 
are given by the same radial wave function, Eq. Q55[l. For pairs with the same I (pseudospin partners) the lower 
components are the same, up to a constant factor. This can be seen from Eqs. (|57|l and l|56l) . since the pseudospin 
partners have quantum numbers (n ,1) (re < 0) and (n— 1,1 + 2) (n > 0). This particular behavior could be related to 
a particular symmetry of the Dirac Hamiltonian in Eq. (J2J) with u>i — 0, as is the case with the pseudospin symmetry, 
to be discussed in Sec. IIVI 



B. Case lo 2 = and A = 



In this subsection we set to zero both A and U potentials. Then, Eq. (|19fl with u>2 = leads to the following second 
order differential equation for the upper component [note that k(k + 1) = 1(1 + 1), irrespective of the sign of n] 



d 2 
d r 2 



1(1 + 1) m(£ + m) 



(to 2 - £ 2 



g K (r) = 0. 



(58) 



The solutions of this equation are given by Eq. ((30(1 . where a is given by Eqs. 1(20(1 - 1(22(1 m which u>2 is set to zero. 
The upper and lower radial functions for two levels corresponding to n with opposite signs are plotted in Fig. EH1 
From this figure we see that the node structure is the same as given in the preceding section when u>2 7^ [Eq. (|34|) ]. 
The energy eigenvalues are readily obtained from Eq. l(2"5)l by just taking the limit w 2 — > 0, yielding 
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Figure 10: Radial wave functions for the potential S(r) = ^mu> 2 r 2 , with [/ = A = for the states 2s]/2, 1^3/2, and 
1(^5/2, respectively, with u>\ — 2 and to = 10. 



2 TO 



( 2n + I + - 



(n = 0,1,2,...). 



(59) 



This equation for the energy, valid for A — 0, shows explicitly that there is no spin-orbit term and that states with 
j = I ± 1/2 are degenerate. We can prove that this equation has only one real solution, and it has to be such that 
£ > to. To begin with, since the right-hand side is always positive, it is clear that any real solution must be greater 
than to, otherwise the left-hand side would be negative. We prove that there is only one of such solutions by squaring 
Eq. (|59|l and using the Descartes' rule of signs. This rule states that the number of positive real roots of an algebraic 
equation with real coefficients a^x k + • • • + a\x + ao = is never greater than the number of changes of signs in the 
sequence a^, ■ ■ ■ , cii, cto (not counting the null coefficients) and, if less, then always by an even number |26ll27|. Since 
we can write the square of Eq. 1)59(1 as 



a 2 = 



(60) 



with x = (£ — m)/m and a = ^/2u\(2n + I + 3/2) /m, then, by Descartes' rule of signs, only one solution with x > 0, 
i.e., £ > to, exists. Since by squaring Eq. 1(59(1 we may only introduce a negative (x < 0) spurious solution, this 
solution is the only positive solution of that equation, and the proof is complete. Therefore, for given values of u)\, n 
and I, i.e, only discrete, positive (greater than to) energies are allowed. 

The nonrelativistic limit is obtained by enforcing u>2 — > in Eq. 1(29(1 . The right-hand side becomes the expression 
for the nonrelativistic energy eigenvalues of the harmonic-oscillator with no spin-orbit coupling. Similarly, Eq. 130(1 
gives the non-relativistic harmonic oscillator radial wave function, with a 2 — tmo\. Figure ITT1 shows the single particle 
spectrum for A = 0, m = 10, and ui\ = 2. For sufficiently small ui\ the energy levels are essentially equidistant, 
a well-known feature of nonrelativistic harmonic-oscillator energy levels, and we do not have spin-orbit interaction: 
states with j = I ± 1/2 are degenerate. These results allow us to stress that the A potential can be closely connected 
to the intrinsic relativistic content of a theory with scalar (S) and vector (V) potentials. The nonrelativistic harmonic 



oscillator is in fact the limit of the relativistic theory with A = V — S = and a potential well E = S + V when uix/m 
becomes very small. 

Thus, an important conclusion of our study is that a Dirac equation with scalar and vector harmonic-oscillator 
potentials with the same sign and magnitude (A = 0), instead of a linear tensor potential U, seems to be the most 
natural way to introduce the harmonic oscillator in relativistic quantum mechanics. 
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Figure 11: Single particle energies for the case A = U = with u)\ = 2, L02 = 0, and m = 10. 



C. Case ^2 = and E = 



Now we consider that both the tensor potential U and £ are zero. The situation is very similar to the preceding 
subsection, except that we now find that it is the lower component that satisfies the following second order differential 
equation 



1(1 + 1) m(8-m) 



2 2 I 2 c2\ 

u)\r — (to — £ ) 



Mr) = 



(61) 



[note that k(k — 1) = k(k + 1) = 1(1 + 1)]. The solution of this equation is the function (|4ri|) setting 0J2 = in Eq. I|37(l . 

These radial wave functions are represented in Fig. El for the 2s\/2, 1^3/2j and 2d 5 / 2 states. Again, we see that 
the node structure is the same as the one given by Eq. (|47l) . For the states 2si/ 2 and ld 3 / 2 , which are pseudospin 
doublets as we will see later, the lower components are equal. This is a general feature of pseudospin symmetry, since 
when one acts upon a given state with the pseudospin SU(2) generators to obtain its pseudospin partner, the radial 
lower component is not changed because of the particular form of these generators, as was shown in Refs. [TfiL Il8j. 

The eigenenergies for S = U = are obtained taking the limit 0J2 — > in Eq. (|41|l . which leads to 



2 TO 



= wi ( 2n + I + - 



(n = 0,1,2,...). 



(62) 



Again, since the right-hand side of Eq. (|62|l is positive and real, we see that the real solutions must have positive 
binding energy E — £ — to. Using similar arguments as with Eq. l|59|l . one can show that there can be just one of 
such solutions. 

The nonrelativistic limit in this case is obtained by setting lo 2 = in Eq. I|42|l . As remarked above, the resulting 
expression is of second order in u>x/m, meaning that the energy is zero up to first order in uji/m. We can interpret 
this fact by saying that, up to this order, there is no nonrelativistic limit for £ = U = and therefore the theory 
is intrinsically relativistic and so is the pseudospin symmetry. In this case, the second order equation in Eq. I|61|l . 
which only depends on I, and also the eigenenergies in Eq. H62|l. show that there is no pseudospin-orbit coupling 
and therefore the states with same (n,l), but with j = I + 1/2 and j = I — 1/2, are degenerate. Thus, when A is 
a harmonic-oscillator potential and S = U = 0, there are only positive-energy bound states and exact pseudospin 
symmetry, i.e., A acts as a binding potential. This is an interesting result in view of the fact that the pseudospin 
symmetry obtained in the limit £(r) — > cannot be realized for nuclear vector and scalar mean fields which go to 
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Figure 12: Radial wave functions for the potential S = of the states (a) lpi/2 (2si/2)> (b) lp 3 / 2 (1^3/2)7 and (c) 
O/5/2 (1^5/2) with the same parameters as in Fig. 1111 

zero as r — ► 00, since in that case E acts as a binding negative central potential well and therefore no bound states 
may exist when E = [T2I Il3j. The spectrum of single particle states for the case E = and m = 10 is shown in 
Fig. llflf a) using the quantum numbers of the upper component, which can be seen as the analog of the nonrelativistic 
quantum numbers. In Fig. ll3f b) we classify the same energy levels by the quantum numbers of the lower components 
f K . The comparison between these two figures manifests the pseudospin symmetry and its quantum numbers. For 
example, the doublets [lsi/2 — 0e? 3 / 2 ] and [lj>3/2 — O/5/2], which have the same pseudo angular momentum I and the 
same n, are, in the new notation, [Opi/2 — OP3/2] and [Oc/3/2 — Odg/a], respectively. Therefore, the harmonic oscillator 
with E = U = and A = ^muifr 2 provides an example of exact pseudospin symmetry. The fact that for this case we 
do not have a nonrelativistic limit, already shows the relativistic nature of the pseudospin symmetry. Figure ITST a) 
shows the (2n + I) degeneracy, which means that not only states with same h, I are degenerate (pseudospin partners) 
but also, for example, (n — 1, 1 + 2) or (n + 1, 1 — 2) have the same energy. 

The singlet states (I — 0) are the states npi/2 with I — 1 inside the small dashed squares in Fig. |21 It is interesting 
to point out that the ground state level for this case is the state Opi/2 with 1 = 1. This is an unusual behavior coming 
from the relativistic nature of this particular case. 

IV. THE PSEUDOSPIN SYMMETRY AND THE INTRUDER LEVELS 

The concept of pseudospin symmetry appeared in nuclear physics more than 20 years ago It was introduced 

to account for the observation, in heavy nuclei, of the quasidegeneracy of orbitals with quantum numbers (for fixed n 
and I) 

(n, l,j = l + ±) and (n - 1, I + 2, j = I + |j . (63) 

Such pairs of single particle states are known as pseudospin partners. 

The doublet structure in the spectrum is better expressed in terms of the pseudo-orbital angular momentum and 
pseudospin quantum number 1 = 1 + 1 and s = s = 1/2, respectively. The former, as noted by Ginocchio |12|. is just 
the orbital angular momentum of the lower component of the Dirac spinor, introduced in Sec. [H] The pseudospin 
partners have the same /. For example, for the partners [nsi/2, (n — 1)^3/2], I — 1, and for [np 3 / 2l (n — 1)^5/2] one 
has I = 2. Yet, the total angular momentum quantum number is still given by j = I ± s, where the aligned (+) and 
unaligned (— ) cases are determined by the sign of R [see Eq. (PU)) ]. Pseudospin partners have k with opposite signs 
(the states with n and n — 1 radial quantum numbers have, respectively, k > and k < 0) and thus pseudospin and 
pseudo-orbital angular momentum have opposite alignments. The fact that I is a good quantum number can then be 
related to the disappearance of the term with E' in Eq. (fTTjl . which can be interpreted as pseudospin-orbit interaction 
term j2lj. From that equation it is also clear that this happens only when {7 = 0, that is, the tensor interaction breaks 
pseudospin symmetry. 
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Figure 13: Single particle energies for the case £ = with the quantum numbers (a) (nlj) and (b) (hi j = j). The 
parameters are lo\ = 2, u>2 = 0, and m = 10. 



The existence of pseudospin partners is connected to a SU(2) symmetry of the Dirac equation with only scalar 
and vector potentials when A = or £ = 0, rather than to the particular shapes of these potentials. Back in 1975, 
Bell and Ruegg 0] obtained the pseudospin generators for that symmetry. In that pioneering work only the case 
A = was studied with the possible application to meson spectra exhibiting a tiny spin-orbit splitting. More recently, 
Ginocchio considered the other case, £ = 0, to explain the quasi-degeneracy of pseudospin doublets in nuclei. The 
generators of the symmetry for radial potentials were worked out in Ref. [16j. 

It is interesting to discuss how the pseudospin symmetry gets broken. The SU(2) generators of the pseudospin 
symmetry are given by 0, 0, El 

S l = Sl -(l-P) + - 2(1+$ = ( a J, (64) 



where 



(T-P (T P 2s p 

c = — —Pi-Si, (65) 



P P P 

and Si = CTi/2 (i — 1, 2, 3). The commutator of this operator with the Hamiltonian in Eq. with U = is 



(66) 



Thus, the pseudospin symmetry breaking can be related to the commutator [£, Sj] |22j. This commutator is zero when 
£ = 0, but for radial potentials it is enough that £'(r) = j2l|. The quasidegeneracy of some pseudospin partners 
can be seen in a Dirac Hamiltonian with scalar and vector potentials of Woods-Saxon type As explained before, 
one cannot set £ = for those kind of potentials, so we cannot get the full degeneracy of all pseudospin partners. In 
fact, the nuclear pseudospin symmetry in nuclei has a dynamical character coming from a cancellation among different 
terms that contribute to the energy as discussed in several works [ljj, |23j, |2J, |2£j . 

We have seen in the previous sections that the relation between the nodes of radial functions depends on whether 
one has A = or £ — in Dirac Hamiltonian with harmonic-oscillator potentials. Leviatan and Ginocchio have 
shown |14| that, in the limit of pseudospin symmetry, the lower radial functions of the pseudospin partners have the 
same number of nodes. Using this information, they were able to show [Lsl ] that, for scalar and vector potentials 
which go to zero at infinity, like Woods-Saxon potentials, the number of nodes is related by 



I "9 > K < 

n f = . , ^ n ( 67 ) 

1 , K > 0. 



This is the behavior we have found before for scalar and vector oscillator harmonic potentials with A = 0, Eq. H34JI . 
precisely the case for which we have the nonrelativistic harmonic-oscillator limit. 

Let us now discuss the levels with n g = nf — 0. From Eq. i|67|) . we see that this can only happen when j — 
I + 1/2 (k < 0) but never with j = I — 1/2 (k > 0). If we recall that the pseudospin doublets have the same n/, in this 
case n/ = 0, states such as Osi/2, OP3/21 0^5/2: etc., have no pseudospin partners. These states are known as intruder 
states. The 0<?g/ 2 state, for example, is an intruder and what would be its pseudospin partner — the state 0in/2 
(k = 6) — does not exist because n g = rif — 1 would be negative for rif = 0. This is what we observe in the nuclear 
spectra. 

In the case of exact pseudospin symmetry (£ = 0) with a harmonic oscillator potential for A, which we have 
discussed in Sec. IIII CI the structure of the nodes it is quite different. In accordance with Eq. 147|) . 

{n„ — 1 , k < 
(68) 
Jig , K > (J. 

Thus the intruder levels that have n g — nt — may only exist in the exact pseudospin limit when j = 1 — 1/2(k > 0), 
as we see from Eq. (I68fl and, in that sense, they are not the intruder levels seen in nuclei, where j = I + 1/2 (k < 0L 

In a recent paper the intruder levels have been discussed in the context of the relativistic harmonic oscillator [Hj 
and it was suggested that they will have partners in the limit of the exact pseudospin symmetry (E = 0). However, 
in that paper, the pseudospin partners have been classified considering only the radial quantum number of the lower 
component (same n/). In that classification, levels with the same rif ~ may have partners according to Eq. I|tj8[) . 
However, it seems that the relation between the number of nodes of the upper and lower components was considered 
to be given by Eq. lj()7(l instead of the correct one given by Eq. 1)6811 . If we take into account that the nuclear intruder 
levels have j = I + 1/2 (k < 0) then, by Eq. (|68[) . we see that states with negative k that have zero nodes in the 
lower component will have a node in the upper component, n g = 1. This disagrees with nuclear spectroscopy, since 
the intruder states should have n g = 0. This is shown in Fig. I13f a) where we show inside the solid square the state 
139/2 with n = 1 and its partner Oin/2 that in the (n, /) classification in Fig. I13f b^ are the pseudospin doublet 
[Ohg/2 - 0h 11/2 \. 

We can summarize this discussion and conclude that even though the harmonic oscillator with E = has an 
exact pseudospin symmetry that could explain the almost degeneracy in I observed in nuclei for the doublets with 
j = I + 1/2 (k < 0), j = I — 1/2 (k > 0), radial and orbital quantum numbers related as shown in Eq. JHSJ, it cannot 
explain the origin of the intruder levels. Furthermore, the harmonic oscillator presents a much higher degeneracy in 
the energy levels besides the simple degeneracy I, as we have shown, which implies that the nuclear potential cannot 
be approximated by a harmonic oscillator and must be more complicated in order to remove this extra degeneracy. 
The asymptotic properties of the central mean-field potential also seem to play a crucial role in defining the node 
structure of the single particle levels. 



V. CONCLUSIONS 



We have presented a generalized harmonic oscillator for spin 1/2 particles that includes not only the usual linear 
tensor potential obtained by a redefinition of the momenta, but also a vector (V) and scalar potential (S) that appear 
in the Dirac Hamiltonian in the specific combinations l/2(l±P)V, in which V is a harmonic oscillator radial potential. 



A special attention was paid to the 1/2(1 — (3)V combination, i.e., when S = — V(H = 0), which has been related 
recently to the nuclear pseudospin symmetry. 

We have derived the energy eigenvalue equations and shown explicitly the wave functions. The nonrelativistic limits 
for all cases have been discussed. The analytical expression for the upper and lower components of the Dirac spinor 
for U = 0, when E = (or A = 0), has been found. 

We have discussed the structure of the solutions of the eigenvalue equations and presented explicitly the positive 
energy solutions. We have concluded that only when the tensor potential is turned on we have the negative bound 
state solutions and in this case, when uj\ = 0, the two bound state solutions are symmetric. When the tensor potential 
is absent we have also shown that the special conditions between the scalar and vector potential (S = V or S — — V) 
needed to have an harmonic-oscillator potential with scalar and vector potentials exclude the negative bound state 
solutions from the spectra. 

The structure of the radial nodes for the Dirac spinor has also been presented and compared to the case when the 
potentials vanish for large distances. We discussed in detail the case of the harmonic oscillator with U = £ = for 
which pseudospin symmetry is exact, i.e., pseudospin doublets (states with same I) are degenerate. In this case we can 
have bound states due to the positive harmonic-oscillator potential A. These states can be regarded as intrinsically 
relativistic, in the sense that the A potential has no analog in the nonrelativistic limit. In fact, they show a very 
peculiar structure of radial nodes in comparison to the nonrelativistic harmonic-oscillator radial wave function: the 
states with j = I + 1/2 have at least one node [n > 1) in the upper component and can have zero nodes only in the 
lower one. As a consequence, the intruder levels on the nuclear spectra (isolated levels with n = and j = I + 1/2) 
do not exist in the harmonic oscillator in the limit of nuclear pseudospin symmetry (X = 0). In this limit the states 
with zero nodes in the upper component have j = I — 1/2. Moreover, the relativistic harmonic oscillator in the case 
U = X = has a much higher energy degeneracy than the approximate pseudospin symmetry seen in the nuclear 
spectra (doublets with same h, I). This is not surprising since the meson exchange theory of nuclear forces produces 
nuclear potentials that vanish when r — > oo, a completely different asymptotic behavior from the harmonic-oscillator 
potential discussed here. One possible way to break this extra degeneracy of the harmonic oscillator is to add a 
Woods-Saxon- like potential as it has been done in Ref. 0] . 

Another important conclusion from our work is that a Dirac equation with scalar and vector harmonic-oscillator 
potentials with the same sign and magnitude (A = 0), instead of the linear tensor potential U, is the most natural 
way to introduce the harmonic oscillator in relativistic quantum mechanics, since one gets the correct nonrelativistic 
limit, so that the energy levels and the upper component of the wave functions are very similar to the non- relativistic 
ones. 

Finally, our analysis of the relativistic harmonic oscillator shows, in a very simple and analytical way, that the 
nuclear pseudospin symmetry (S = 0) does not have a nonrelativistic limit, in the sense that the eigenvalues are of 
second order in a Lo/m expansion. So we conclude that the relativistic harmonic oscillator is an excellent example 
to show explicitly what Ginocchio had already emphasized in his pioneering work: the pseudospin symmetry in the 
Dirac Hamiltonian (when S = 0) is indeed a symmetry of relativistic nature. 
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